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Abstract 

We study the complexity of approximately solving the weighted counting constraint 
satisfaction problem #CSP(J r ). In the conservative case, where J- contains all unary 
functions, there is a classification known for the case in which the domain of functions 
in T is Boolean. In this paper, we give a classification for the more general problem 
where functions in J- have an arbitrary finite domain. We define the notions of weak 
log-modularity and weak log-supermodularity. We show that if T is weakly log-modular, 
then #CSP(J r ) is in FP. Otherwise, it is at least as difficult to approximate as #BIS, 
the problem of counting independent sets in bipartite graphs. #BIS is complete with 
respect to approximation-preserving reductions for a logically defined complexity class 
^RHIIi, and is believed to be intractable. We further sub-divide the #BIS-hard case. If 
T is weakly log-supcrmodular, then we show that #CSP(J 7 ) is as easy as a Boolean log- 
supermodular weighted #CSP. Otherwise, we show that it is NP-hard to approximate. 
Finally, we give a full trichotomy for the arity-2 case, where #CSP(J r ) is in FP, or is 
^BIS-cquivalent, or is equivalent in difficulty to #SAT, the problem of approximately 
counting the satisfying assignments of a Boolean formula in conjunctive normal form. We 
also discuss the algorithmic aspects of our classification. 



1 Introduction 

In the weighted counting constraint satisfaction problem, there is a fixed finite domain D 
and a fixed finite "weighted constraint language" J-, which is a set of functions. Every 
function F £ T maps a tuple of domain elements to a value called a "weight". In the 
computational problem #CSP(J 7 ), an instance consists of a set V = {vi, . . . , v n } of variables 
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and a set of "weighted constraints" . Each weighted constraint applies a function from F to 
an appropriate-sized tuple of variables. 

For example, with the Boolean domain D = {0, 1} we could consider the situation in which 
F consists of the single binary (arity-2) function F defined by F(0, 0) = F(0, 1) = F(l, 0) = 1 
and F(l, 1) = 2. We can construct an instance with variables v±, V2 and v% and weighted 
constraints F{v\,V2) and F(v2,v^). If x is an assignment of domain elements to the variables 
then the total weight associated with x is the product of the weighted constraints, evaluated 
at x. 

For example, the assignment that maps v±, V2 and v$ all to has weight F(0, 0)F(0, 0) = 1 
but the assignment that maps all of them to 1 has weight F(l, 1)F(1, 1) = 4. Two assignments 
have weight F(l, 1).F(1, 0) = -F(0, 1).F(1, 1) = 2 and the other four assignments each have 
weight 1. The computational problem is to evaluate the sum of the weights of the assignments. 
For this instance, the solution is 13. 

There has been a lot of work on classifying the computational difficulty of exactly solving 
#CSP(J r ). For some weighted constraint languages F, this is a computationally easy task, 
while for others, it is intractable. We will give a brief summary of what is known. For more 
details, see the surveys of Chen [12J and Lu |26j. 

First, suppose that the domain D is Boolean (that is, suppose that D = {0,1}). For 
this case, Creignou and Hermann [15] gave a dichotomy for the case in which weights are 
also in {0, 1}. In this case, they showed that #CSP(F) is in FP (the set of polynomial-time 
computable function problems) if all of the functions in F are affine, and that otherwise, it is 
#P-complete. Dyer, Goldberg, and Jerrum [T7] extended this to the case in which weights are 
non-negative rationals. For this case, they showed that the problem is solvable in polynomial 
time if (1) every function in F is expressible as a product of unary functions, equalities and 
disequalities, or (2) every function in J 7 is a constant multiple of an affine function. Otherwise, 
they showed the problem is complete in the complexity class FP* P . We will not deal with 
negative weights in this paper. However, it is worth mentioning that these results have been 
extended to the case in which weights can be negative [5], to the case in which they can be 
complex [23], and to the related class of Holant* problems [11]. Other dichotomies are also 
known for Holant problems (see [26J). 

Next, consider an arbitrary finite domain D. For the case in which weights are in {0, 1}, 
Bulatov's breakthrough result [3] showed that #CSP(J-") is always either in FP or #P-hard. 
A simplified version was given by Dyer and Richerby [19|,l20j who introduced a new criterion 
called "strong balance" . The dichotomy was extended to include non-negative rational weights 
by Bulatov et al. [I] and then to include all non-negative algebraic weights by Cai, Chen and 
Lu [9l[T0]. Cai, Chen and Lu gave a generalised notion of balance that we will use in this 
work. Finally, Cai and Chen [8] extended the dichotomy to include all algebraic complex 
weights. The criterion for the unweighted #CSP dichotomy is known to be decidable [20J and 
this carries through to non- negative rational weights and non-negative algebraic weights [9]. 
Decidability is currently open for the complex case. 

Much less is known about the complexity of approximately solving ^CSP(J-"). Before 
describing what is known, it helps to say a little about the complexity of approximate counting 
within ^tP. Dyer, Goldberg, Greenhill and Jerrum [16] identified three complexity classes for 
approximation problems within #P. These are: 

1. problems that have a fully polynomial randomised approximation scheme (FPRAS), 

2. a logically defined complexity class called #RfflIi, and 
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3. a class of problems for which approximation is NP-hard. 

A typical complete problem in the class ^RHQi is #BIS, the problem of approximately 
counting independent sets in bipartite graphs. It is known that either all complete problems 
in ^RHIIi have an FPRAS, or none do; it is conjectured that none do [21|.I22| . A typical 
complete problem in the third class is ^SAT, the problem of counting satisfying assignments 
of a Boolean formula in conjunctive normal form. Another concept that turns out to be 
important in the classification of approximate counting CSPs is log-supermodularity [HE]. A 
function with Boolean domain is log-supermodular if its logarithm is supermodular; we give 
a formal definition later. 

Given those rough ideas, we can now describe what is known about the complexity of 
approximately solving #CSP(J r ) when the domain, D, is Boolean. For the case in which 
weights are in {0, 1}, Dyer, Goldberg and Jerrum gave a trichotomy [18| . If every function in 
J- is affine, then #CSP(J~) is in FP. Otherwise, it is as hard to approximate as #BIS. The 
hard approximation problems are divided into #BIS-equivalent cases (which arise when the 
functions in J- can be expressed using "implies" and fixing the values of certain variables) 
and the remaining cases, which are all shown to be # SAT-equivalent. 

In the more general case where the domain D is still Boolean, but the weights can be 
arbitrary non-negative values, no complete classification is known. However, Bulatov et al. [7J 
were able to give a classification for the so-called "conservative" case, in which T contains 
all unary functions. Their result is reproduced as Lemma [7J below. Here is an informal 
description. 

• If every function in T can be expressed in a certain simple way using disequality and 
unary functions, then, for any finite Q C T, #CSP(C/) has an FPRAS. 

• Otherwise, 

— there is a finite Q C J- such that #CSP(£7) is at least as hard to approximate as 
#BIS and, 

— if J- contains any function that is not log-supermodular, then there is a finite 
Q C T such that #CSP(£/) is at least as hard to approximate as #SAT. 

Yamakami [32] has also given an approximation dichotomy for the case in which even more 
unary functions (including those with negative weights) are assumed to be part of J 7 . The 
negative weights introduce cancellation, making more weighted constraint languages J- in- 
tractable. In this paper, we stick to the non-negative case, in which more subtle complexity 
classifications arise. 

Prior to this paper, there were no known complexity classifications for approximately 
solving #CSP(J-") for the case in which the domain D is not Boolean. Thus, this is the 
problem that we address in this paper. Our main result (Theorem below) is a complexity 
classification for the conservative case (where all unary weights are contained in J 7 ). Here is 
an informal description of the result. 

• If J 7 is "weakly log-modular" (a concept we define below) then, for any finite Q C J-, 
#CSP(£) is in FP. 

• Otherwise, there is a finite Q C T such that ^CSP(^) is at least as hard to approximate 
as #BIS. Furthermore, 
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— if J 7 is "weakly log-supermodular" (again, defined below) then, for any finite Q C F, 
there is a finite set Q' of log-supermodular functions on the Boolean domain such 
that #CSP(£/) is as easy to approximate as #CSP(£/'); 

— otherwise, there is a finite Q C F such that #CSP(£/) is as hard to approximate 
as #SAT. 

Informally, F is weakly log-supermodular if, for every binary function F that can be 
expressed using functions in F, every projection of F onto two domain elements is log- 
supermodular (see Definition Thus, in some sense, our result shows that all the diffi- 
culty of approximating conservative weighted constraint satisfaction problems arises in the 
Boolean case. Even when the domain D is larger, approximations which are #SAT-equivalent 
are #SAT-equivalent precisely because of intractable Boolean problems which arise as sub- 
problems. 

In addition to the complexity classifications described above (FP versus #BIS-hard and 
"as easy as a Boolean log-supermodular problem" versus #SAT-equivalent) we also give a 
full trichotomy for the binary case (i.e., where all functions in T have arity 1 or 2). 

• If J 7 is weakly log-modular then, for any finite Q C F, #CSP(£?) is in FP. 

• Otherwise, if F is weakly log-supermodular, then 

— for every finite Q C F, #CSP(£?) is as easy to approximate as #BIS and 

— there is a finite Q C F such that #CSP(<5) is as hard to approximate as #BIS. 

• Otherwise, there is a finite Q C F such that #CSP((/) is as hard to approximate as 
#SAT. 

The final section of our paper discusses the algorithmic aspects of our classification for 
the case in which F is the union of a finite, weighted constraint language Ti and the set of 
all unary functions. In particular, we give an algorithm that takes T-L as input and correctly 
makes one of the following deductions: 

1. #CSP(£) is in FP for every finite G C F; 

2. #CSP(Q) is LSM-easy for every finite Q C F and ^BIS-hard for some such Q; 

3. #CSP(£?) is #BIS-easy for every finite Q C F and #BIS-equivalent for some such Q; 

4. #CSP(<5) is #SAT-easy for all finite Q C F and # SAT-equivalent for some such Q. 
Further, if every function in % has arity at most 2, the output is not deduction 

1.1 Previous work 

The first contribution of our paper is to show that, if F is weakly log- modular then, for any 
finite Q C F, #CSP(£) is in FP. Otherwise, there is a finite Q C F for which #CSP(£) 
is at least as hard to approximate as #BIS. We also show that, if F is not weakly log- 
supermodular, then there is a finite Q C F, such that #CSP(£) is #SAT-hard. This work is 
presented in Sections [2] and [3] below and builds on two strands of previous work. 
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• The hardness results build on the approximation classification in the Boolean case [7J 
and, in particular, on the key role played by log-supermodular functions. 

• The easiness results build on the classification of the exact evaluation of #CSP(J-") in 
the general case [9], and in particular on the key role played by "balance". 

The second contribution of our paper is to show that, if J- is weakly log-supermodular, 
then, for any finite Q C J 7 , there is a finite set Q' of log-supermodular functions on the Boolean 
domain such that #CSP(£/) is as easy to approximate as #CSP((7'). This builds on three 
key studies of the complexity of optimisation CSPs by Takhanov [29.30J, Cohen, Cooper and 
Jeavons [13] and Kolmogorov and Zivny |24|125] . In all three cases, we use their arguments 
and ideas, and not merely their results. Thus, we delve into these three papers in some detail. 

Our final contribution is the trichotomy for the binary case. This relies additionally on 
work of Rudolf and Woeginger [28] on decomposing matrices known as Monge matrices. 

1.2 Preliminaries and statement of results 

Let D be a finite domain with \D\ > 2. It will be convenient to refer to the set Funcfc(L>,i?) 
of all functions D k — > R for some codomain R, and the set Func(.D, R) = UfcLo Funcfc(D, R). 
Let EQ be the binary equality function defined by EQ(x, x) = 1 and EQ(x, y) = for x ^ y; 
let NEQ(x,y) = l-EQ(x,y). 

We use the following definitions from [JJ. Let J 7 be a subset of Func(D,i?). Let V = 
{vi, . . . , v n } be a set of variables. An atomic formula has the form tp = G{x>i 1 , . . . , Vi a ) where 
G € J 7 , a = a(G) is the arity of G, and [vi x ,Vi 2 , . . . , vt a ) € V a is called a "scope". Note that 
repeated variables are allowed. The function F v : D n — > R represented by the atomic formula 
(p = G(u iu . . . ,v ia ) is just ^(x) = G(x(wj 1 ), . . . ,x(v ia )), where x: {vi, . . . , v n } -s> D is an 
assignment to the variables. To simplify the notation, we write Xj = x(i>j) so 

-?^(x) = Gixi^ %i a )- 

A pps-formula ( "primitive product summation formula" ) is a finite summation of a finite 
product of atomic formulas. A pps-formula ip over T in variables V = {v\, . . . , v n+ k\ has the 
form 

m 

V n + l,...,V n + k j = l 

where (pj are all atomic formulas over T in the variables V. (The variables V are free, and the 
others, V' \ V, are bound.) The formula tp specifies a function : D n — > R in the following 
way: 

m 

yeD k j=i 

where x and y are assignments x: {vi, . . . ,v n } — > D and y : {v n+ i, . . . ,v n+ k} — > D. The 
functional clone (J~)# generated by T is the set of all functions that can be represented by a 
pps-formula over T U {EQ}. Crucially, ((~F)#)# = {J~)# [3 Lemma 1]; we will rely on this 
transitivity property implicitly. 

Definition 1. A weighted constraint language J 7 is a subset of Func(D, Q>o). Functions in 
J- are called weight functions. 
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In Section|3]we will introduce valued constraint languages and cost functions, which pertain 
to optimisation CSPs. It is important to distinguish these from the weighted version, which 
is used for counting. 

Definition 2. A weighted constraint language T is conservative if Ud C F, where Uo = 
Funci(D,Q> ). 

Definition 3. A weighted constraint language J- is weakly log-modular if, for all binary 
functions F € {P)# and elements a, b € D, 

F(a,a)F(b,b) = F(a,b)F(b,a), or 
F(a,a) = F(b,b) = 0, or 

F(a, b) = F(b, a) = 0. (1) 

Definition 4. T is weakly log-supermodular if, for all binary functions F £ (F)# and elements 
a,b £ D, 

F{a,a)F(b,b) > F(a,b)F(b,a) or F(a, a) = F(b, b) = 0. (2) 

Definition 5. A function i* 1 € Funcfc({0, l},Q>o) is log-supermodular if 

F(xVy)F(xAy) > F(x)F(y) 

for all x, y E {0, l} k , where A (min) and V (max) are applied component- wise. LSM is the set 
of all log-supermodular functions in Func({0, l},Q>o). 

It is known Lemma 7] that (LSM)^t = LSM. Here is a precise statement of the computa- 
tional task that we study. A #CSP problem is parameterised by a finite, weighted constraint 
language T as follows. 

Name #CSP(J"). 

Instance A pps-formula if) consisting of a product of m atomic F-formulas over n free vari- 
ables x. (Thus, ij) has no bound variables.) 

Output The value ^2 x€D n F/,(x) where F^ is the function defined by ifi. 

Where convenient, we abuse notation by writing ^CSP(-F) to mean #CSP({-F}) and by 
writing #CSP(F, P) to mean #CSP(TUP). 

As in [7] (and other works) we take the size of a #CSP(J-") instance to be n + m, where 
n is the number of (free) variables and m is the number of weighted constraints (atomic 
formulas). In unweighted versions of CSP and #CSP, we can just use n as the size of an 
instance, since the number of constraints can be bounded by a polynomial in the number of 
variables. However, in weighted cases, the multiplicity of constraints matters so we cannot 
bound m in terms of n. We typically denote an instance of #CSP(J-") by / and the output 
by Z(I), which is often called the "partition function". 

A counting problem, for our purposes, is any function from instances (encoded as words 
over a finite alphabet X) to Q>o- A randomised approximation scheme for a counting prob- 
lem #X is a randomised algorithm that takes an instance w and returns an approximation 
Y to #X{w). The approximation scheme has a parameter e G (0, 1) which specifies the error 
tolerance. Since the algorithm is randomised, the output Y is a random variable depending 
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on the "coin tosses" made by the algorithm. We require that, for every instance w and every 
£6(0,1), 

Pr [e- £ #X(w) < Y < e £ #X(w)] > 3/4. (3) 

The randomised approximation scheme is said to be a fully polynomial randomised approxi- 
mation scheme, or FPRAS, if it runs in time bounded by a polynomial in \w\ (the length of 
the word w) and See Mitzenmacher and Upfal \27\ Definition 10.2]. 

Suppose that f^X and #Y are two counting problems. An "approximation-preserving 
reduction" (AP-reduction) [16J from j^X to #Y gives a way to turn an FPRAS for j^Y 
into an FPRAS for Specifically, an AP-reduction from f^X to f^Y is a randomised 

algorithm A for computing ^X using an oracle for #Y. The algorithm A takes as input 
a pair (w, e) S X* x (0,1), and satisfies the following three conditions: (i) every oracle call 
made by A is of the form (v, 5), where v £ X* is an instance of #Y, and < S < 1 is an error 
bound satisfying d^ 1 < poly(|w|, e^ 1 ); (ii) the algorithm A meets the specification for being 
a randomised approximation scheme for j^X (as described above) whenever the oracle meets 
the specification for being a randomised approximation scheme for j^Y; and (hi) the run- 
time of A is polynomial in \w\ and e^ 1 . If an AP-reduction from #X to j^Y exists we write 
#X <ap jfY . Note that, subsequent to [16], the notation <ap has been used to denote a 
different type of approximation-preserving reduction which applies to optimisation problems. 
In this paper, our emphasis is on counting problems so we hope this will not cause confusion. 

The notion of pps-definability described earlier is closely related to AP-reductions. In 
particular, Lemma 17] shows that G € implies that #CSP(J r ,G) < A p #CSP(J"). 

We will use this fact without comment. 

As mentioned above, #BIS is the problem of counting the independent sets in a bipartite 
graph and #SAT is the problem of counting the solutions to a Boolean formula in conjunctive 
normal form. We say that a counting problem j^X is #Y-easy if j^X <ap ifY and that it is 
#Y-hard if j^Y <ap f^X. A problem j^X is LSM-easy if there is a finite, weighted constraint 
language T C LSM such that #A < AP #CSP(J r ). 

WE now state our main theorem. Note that we have only defined the problem #CSP(J-") 
for finite languages whereas conservative languages are, by definition, infinite. 

Theorem 6. Let T be a conservative weighted constraint language taking values in Q>o- 

1. If J- is weakly log-modular then #CSP(£?) is in FP for every finite Q C J- '. 

2. If J- is weakly log-supermodular but not weakly log-modular, then #CSP(G) is LSM-easy 
for every finite Q C T and #BIS-hard for some such Q . 

3. If J- is weakly log-supermodular but not weakly log-modular and consists of functions 
of arity at most two, then #CSP(<5) is #BIS-easy for every finite Q C J- and #BIS- 
equivalent for some such Q . 

4- If T is not weakly log-supermodular, then #CSP(£) is ^SAT-easy for every finite Q C T 
and #S AT -equivalent for some such Q . 

In particular, among conservative #CSPs, there are no new complexity classes below 
#BIS or above LSM; furthermore there is a trichotomy for conservative weighted constraint 
languages with no functions of arity greater than two. 

The $T3IS-hardness and #SAT-equivalence are proved in Section [2] where they are re- 
stated as Theorem [TUJ The membership in FP is established as Theorem [15] at the end of 
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Section [3J LSM-easiness and $L3IS-easiness are established by Theorem H7J at the end of 
Section [BJ Algorithmic aspects are discussed in Section [TJ 

2 Hardness results 

Our hardness results use the following result from [7J. 

Lemma 7. P7L Theorem 18] Let T be a finite, weighted constraint language with D = {0, 1}. 

• IfTc (NEQ,U{o,i})# then, for any finite S C U {0A} , #CSP(J r U S) has an FPRAS. 

• If T <£. (NEQ,W{o,i})#, then there is a finite S C U{o.i} such that ^CSPp 7 U S) is 
#BIS-hard. 

• If J- <f- (NEQ,W{o,i})# and T (£ LSM, then there is a finite S C ^{0,1} such that 
#CSP(FUS) is#SAT-hard. 

Remark 8. In the statement of [7, Theorem 18], the set Z//{o,i} is replaced with B\, the set 
of all unary functions from {0, 1} to the set of non-negative efficiently computable reals. In 
this paper, we restrict to rationals for simplicity. Even though the statement of [7J Theorem 
18] does not imply Lemma [7J the proof of Theorem 18] does establish the lemma. No 
functions in B\ with irrational weights are used explicitly in the proof — unary functions 
that are used (for example, in the proof of Lemma 12]) are constructed by multiplying 
and dividing other values in the codomains of functions in T . 

In fact we will only use the following special case of Lemma [JJ 

Lemma 9. P7| Theorem 18] Let F be a function in Func2({0, l},Q>o). 

• IfFi (NEQ,W {0il} ) # then {F} UW {0il} is #BlS-hard. 

• IfFi (NEQ,^ {0)1} ) # U LSM then {F}UU {0A} is #SAT-hard. 
Our hardness results now follow from Lemma 

Theorem 10. Let J 7 be a conservative weighted constraint language taking values in Q>o- 

• If J- is not weakly log-modular, there is a finite Q C T such that #CSP(<5) is #BIS-hard. 

• If F is not weakly log-supermodular, there is a finite Q C T such that #CSP(£7) is 
#SAT-hard. 

• For all finite Q C T, #CSP{Q) is #SAT-easy. 

Proof. First, we establish the hardness results. 

Suppose that T is not weakly log- modular. Let H E (J r }# be a function violating (TTJ 
and let a and b be the relevant domain elements, which must be distinct. Let ip: {0, 1} — > D 
be a unary function with c/?(0) = a and ip(l) = b. Define H^: {0, l} 2 — > Q>o by H v (x,y) = 
H((p(x),(f(y)). The following three equations must all fail to hold: 

^(0,0)^(1,1) = ^(0,1)^(1,0) 
^(0,0) = ^(1,1) =0 
H ip (0,l) = H v (l,0) = 0. 
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By Remark 14], every binary function in (NEQ, W{ 0j i})# has one of three forms: Ui(x)U2(y), 
U(x)EQ(x,y) or U (x)NEQ(x , y) . Therefore, H v <£ (NEQ, W {0 ,1}}#- % Lemma there is a 
finite set S C W {0 ,i} such that # BIS ^ap #CSP(^, 5). 
For each U G Z^{o,i}! define t/^-i G kfo by 



tL,-i(x) 



C/(0) if x = a, 
U(l) ifx = b, 
otherwise. 



Let F(0) = F(l) = 1. Let 5' = {U v -i \ U G S U {F}}. Note that {#} U 5' C {F,U D )# is 
finite. 

We describe a reduction from #CSP(/f ¥ ,,5) to #CSP(iJ, 5"). Given an instance I of 
#CSP(H,p, S), replace each use of F^ by F, and each use of C/ G 5 by U^-i G 5", and 
introduce an atomic formula EL-i('w) for each variable to obtain a new instance /' of 
#CSP(F,S') with Z(J) = Z(F). Thus #CSP(F,S / ) is #BIS-hard. 

A similar argument shows that T is #SAT-hard if it is not weakly log-supermodular. In 
this case, we start with a function H G (J-)# violating ([2]) on the elements a,b G D. Defining 
ip and as above, we find that G" LSM. Since H also violates ([I]) on a, b, the argument 
above establishes H v ^ (NEQ, W{o,i})#- By Lemma[9]there is a finite set S C ^{0,1} such that 
#SAT < AP #CSP(^, 5). We then proceed as before. 

# SAT-easiness follows from the construction in Section 3 of [16j . which shows that any 
problem in #P is #SAT-easy. The weighted counting CSPs that we deal with here are 
equivalent, by [4j, to unweighted ones, which are in #P. □ 



3 Balance and weak log-modularity 

In this section we show that weak log-modularity implies tractability, by showing that every 
weakly log-modular weighted constraint language is balanced in the following sense. 

We may associate a matrix M with an undirected bipartite graph Gm whose vertex 
partition consists of the set of rows R and columns C of M. A pair (r, c) G R x C is an edge 
of Gm if) and only if, M rc 7^ 0. A block of M is a submatrix whose rows and columns form 
a connected component in Gm- M has block-rank 1 if all its blocks have rank 1. 

We say that a weighted constraint language T is balanced [§] if, for every function 
F(xi, . . . ,x n ) G (F)# with arity n > 2, and every k with < k < n, the \D\ k x \D\ n ~ k 
matrix F((xi, . . . , (x^+i, . . . ,x n )) has block-rank 1. (This notion reduces to Dyer and 
Richerby's notion of "strong balance" [20] in the unweighted case.) 

A function F : {0, l} n — > R has rank 1 if it has the form F(x±, . . . , x^) = U\{x\) ■ ■ ■ Uk{xk)- 
Given a non-singular two- by- two matrix T G R 2x2 (possibly negative- valued) we define 
T® n F: {0,l} n ->R by 

(T® n F)(xi, ...,x n )= ^ (f[T Xiy ) F(y u ...,y n ) 
y lt ...,y n e{0,l} \i=l / 

The rows and columns of T are considered to be indexed by {0, 1}. We associate with any 
function F : {0, l} 2 ->■ R, the matrix M F G R 2x2 defined by (Afpjy = F(i,j). 

Lemma 11. lei M G R fcxfe . lei F: {0, l} k R. Zei T G R>q 2 6e non-singular. 



9 



1. 



A 2 x 2 matrix M has block-rank 1 if and only if it has rank 1 or it has at most two 
non-zero entries. F: {0, l} 2 — > M has rank 1 if and only if detMp = 0. 



2. M has block-rank 1 if and only if the matrix 

%,u,u'v,v' ; 



M(u,v) M(u,v') 
M(u',v) M(u',v') 



has block-rank 1 for every u, u', v, v'. 



3. (Topkis's theorem) If F is strictly positive and not of rank 1, there is a function 
F' : {0, l} 2 — > R of the following form that is not of rank 1. 

F (xj , Xj) — F(c\ , . . . , Cj_i , Xj, Cj+i , . . . , Cj—\ j Xj , Cj-f-i ; • • ■ j Cfc). 

Here 1 < i < j < k, and each q is a fixed element of {0, 1}. 

4- F has rank 1 if and only ifT® k F has rank 1. 

Proof. 1. A 2 x 2 matrix that has block-rank 1 either has rank 1 or is diagonal or anti- 
diagonal so has two zeroes. Conversely, a matrix that has rank 1 has no submatrix 
whose rank exceeds 1, so has block-rank 1. A matrix with two or more zeroes has no 
2x2 block so can only have blocks of rank 1. 

For the second statement, if F has rank 1 then there are unary functions Uq and U\ so 
that F(x,y) = Uo(x)Ui(y), which implies that detMp = 0. Going the other way, if F 
is identically then it has rank 1. Otherwise, suppose F(i,j) ^ 0. Let Uq{x) = F(x,j) 
and Ui(y) = F(i,y)/ F(i, j). If detM F = then F(x,y) = Uo{x)U x (y), so F has rank 1. 

2. [201 Lemma 38]. 

3. Say that a strictly positive function F is log-modular if / = log F is modular: that is, 
-F(x V y)i ? (x A y) = i ? (x)F(y) for all x, y € {0, l} k . A modular function is an affine 
map (see for example [2j Proposition 24]), so a strictly positive log-modular function is 
a product of unary functions, so it has rank 1. The result is then Topkis's theorem [31] 
in the form stated in [7., Lemma 8]. 

4. If F is of the form Ui(xi) ■ ■ ■ U n {x n ) then 

(T®»F)(zi, • • • , O = (T^C/iXxi) • • • (T m U n )(x n ) 

The reverse implication follows from (^T^ 1 )^ n T^ n F = F, where T~ l is the matrix 
inverse of T. 

□ 

A function F: D n — > Q>o is essentially pseudo-Boolean if its support is contained in a 
set D\ x • • • x D n with \D\\, . . . , \D n \ < 2. The projection of a relation R C D n onto indices 
1 < i < j < n is the set of pairs (a, b) G D 2 such that there exists x € R with xi = a and 
Xj = b. A generalised NEQ is a relation of the form {(xi,Xj), (yi,yj)} C D 2 for some Xi ^ yi 
and Xj ^ yj. 
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Lemma 12. Let F: D n — > Q>o be an essentially pseudo-Boolean function which is not of 
rank 1, and assume that no binary projection of the support of F is a generalised NEQ. Then 
{F} U Ud is not weakly log-modular. 

Proof. Let the support of F be contained in D\ x • • • x D n where \D{\ = 2 for all i. Choose 
bijections pi : {0, 1} — > A for each 1 < i < n. Define F p : {0, l} n — > Q>o by 

F p (x 1 , . . . ,x n ) = F(p 1 (x 1 ), . . .,p n {x n )) 

for all x%, . . . ,x n G {0, 1}. Let T = {\\) an d note that T® n F p is strictly positive. Since F 
is not of rank 1, F p is not of rank 1, so by Lemma [TT1 part Q, T® n F p is not of rank 1. By 
Lemma [TT1 part ([3]), there is a function B: {0, l} 2 — > Q>o of the following form that is not of 
rank 1. 

B(xi, ■Ej ) — -^p) (ci j • • • j Cj— 1 j 3<i j Cj+1 > • • • j Cj— 1) j Cj+i , . . . , C n ) . 

For all indices k G {1, . . . , n} \ {i,j}, define U k G by U k {pk{xk)) = T CkXk for all x fc G 
{0,1}, and ?7 fc (2) = if z £ D k . Define G,H:D 2 ^ Q> and G pi;Pj , H pitPi : {0, l} 2 -> Q> as 
follows. Note in these definitions that i and j are fixed, but pi and are used as subscripts 
in the name of some of the functions as a reminder of the bijections that are being applied 
to the inputs. Thus, in H pi)Pi , the bijection pi is applied to both arguments, even though the 
function depends on both i and j. 

for all yi,yj G D 

for all Xi, Xj G {0, 1} 
for all y', y" G D 
for all x',x" G {0,1} 

where the first sum is over all y\, . . . , Jft+i, • • • , yj-\,yj+\, . ■ ■ ,y n G .D and the second 
sum is over all xi, . . . , . . . , . . . ,x n G {0, 1}. 

Note that M Hpj . Pi = M Gpi P ,M*, = T^M^r^^T^M^T- 1 )* where t denotes 

transpose. Taking determinants and applying Lemma [11] part ([1]) this implies that H puPi 
is not of rank 1. Also, since T is strictly positive, the support of G Pi)Pj is the binary pro- 
jection of the support of F p onto i and j which, by assumption, is not NEQ or EQ. Hence 
H piPi is strictly positive but not of rank 1. Again using Lemma [TT1 part ([TJ we see that H is 
a witness that {F} U Up is not weakly log-modular. □ 

Lemma 13. Every conservative weakly log-modular weighted constraint language is balanced. 

Proof. Let J 7 be a conservative weighted constraint language that is not balanced. We will 
show that T is not weakly log- modular. 

By the definition of balance, there is a function F G {F)# of arity n and a partition 
x = (u, v) of its n variables, such that the matrix F(u, v) is not of block-rank 1. By LemmafTTI 
part ([2]) there is a two- by- two submatrix = Np u u / v v / that is not of block-rank 1. 



i,---,y n ) 



Gpi,pj{xi,Xj) — I T Ck ,x k J Fpipli ■ ■ ■ j x n) 



H(y',y") = ^G(y',y)G(y",y) 



y&D 



H Pi , Pi {x\x") — Gp u p j (x',x)Gp ijPj (x",x) 
xe{o,n 
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Construct an essentially pseudo-Boolean function G from F as follows. For all 1 < i < n 
let Ui G (Md)# C (F)# be the indicator function of D 1 ^ 1 x D f x L> n ~*, where A = {uj, u^} for 
all 1 < i < k and Di = {vi^v[} for all k < i < n. Let G = F\\^ =l Ui. Then NQ u u t v y = iV 
is not of block-rank 1, and G is essentially pseudo-Boolean. 

If the binary projection of the support of G onto two indices i,j is a generalised NEQ 
{(xi, Xj), (yi, Uj)}, construct (G', p(u), p(u'), p(v), p(v')) from (G, u, u', v, v') as follows. Let 
p: D n — > D 11 ^ 1 be the projection operator sending x to x\, . . . , Xj_i, Zj+i, • • • , x„ and let 

G '( x ) = E P (x')=x G ( x ') for a11 x G Note that ' for a11 x G Dn > G ( x ) # G '(P( X )) 

implies that G(x) = because G(x) = unless Xi ^ Xj. Note that N has at least 

three non-zero entries by Lemma [11] part ([T|). So the corresponding three pairs out of 
((u,v)i,(u,v)j), ((u,v')i, (u,v')j), ((u' v)i, (u', v)j), and ((u' v')j, (u' v') 3 -) must each be ei- 
ther (x^ Xj) or (yi,yj). But then the fourth pair must also be (xj, Xj) or (yi,Uj), which implies 
that Nqi p r u \ pf u /\ p r y \ pf^/s = N. Also, G' is essentially pseudo-Boolean, and G' is obtained by 
summing the i'th variable, so G' G (G)#. 

Repeating this process if necessary, we obtain (G', x, x', y, y') such that G' is an essentially 
pseudo-Boolean function in {F,Ud)# = (•7 7 )# and none of the binary projections of the 
support of G' is a generalised NEQ, and Ng',*. n',y,y' is n °t of block-rank 1. So, in particular, 
G' is not of rank 1. By Lemma fT2l {G'} UUd is not weakly log-modular, so (F)# is not 
weakly log-modular. □ 

We now return to Theorem [6] and prove the tractable case. The proof relies on an impor- 
tant theorem of Cai, Chen and Lu about the complexity of exact evaluation. 

Lemma 14. Let F be a finite, weighted constraint language taking non-negative algebraic 
real values. If F is balanced, then #CSP(F) is in FP, and otherwise #CSP(F) is jfP-hard. 

Theorem 15. Let F be a conservative weighted constraint language taking values in Q>o- If 
F is weakly log-modular then, for any finite Q C F, #CSP(ty) € FP. 

Proof. By Lemma [T3l F is balanced. Hence, every finite Q C F is balanced, which implies 



4 Valued clones, valued CSPs and relational clones 

To define valued clones, we use the same set-up as Section 11.21 except that summation is 
replaced by minimisation and product is replaced by sum. Let D be a finite domain with 
\D\ > 2 and let R be a codomain with {0, oo} C R y where oo obeys the following rules for all 
xG-R: x + oo = oo, x<oo and min{x, oo} = x. Let <I> be a subset of Func(.D, R) and let 
V = {vi, . . . , v n } be a set of variables. For each atomic formula 92 = Giy^ , . . . , Vi a ) we use 
the notation f v to denote the function represented by ip, so /^(x) = G(xi 1 , . . . ,Xj Q ). 

A psm- formula ("primitive sum minimisation formula") is a minimisation of a sum of 
atomic formulas. A psm-formula ij) over $ in variables V = {vi, . . . , v n+ k} has the form 



that #CSP(£) is in FP by Lemma QU 



□ 



m 




(4) 
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where tpj are all atomic formulas over $ in the variables V. The formula ip specifies a function 
: D n — > R in the following way: 



(x) = min V f v (x, y) , (5) 

y eD j= i 



where x and y are assignments x: {v\, . . . ,v n } — > D and y : {v n+ \, . . . ,i> n +fc} — ► D. 

The valued clone (&)v generated by $ is the set of all functions that can be represented 
by a psm- formula over <& U {eq}, where eq is the binary equality function on D given by 
eq(x, x) = and eq(a;, y) = oo for x y. 

We next introduce valued constraint satisfaction problems (VCSPs), which are optimisa- 
tion problems. In the work of Kolmogorov and Zivny [24J, the codomain is R = Q>o U {oo}. 
For reasons which will be clear below, it is useful for us to extend the codomain to include 
irrational numbers. This will not cause problems because, with the exception of Theorem [371 
we use only formal calculations from their papers, not complexity results. For Theorem l37l we 
avoid irrational numbers and, in fact, restrict to cost functions taking values in {0, oo} C R. 
Furthermore, all the real numbers we use are either rationals or their logarithms so are effi- 
ciently computable. 

Let M>o = M>o U {oo} be the set of non-negative real numbers together with oo. 

Definition 16. A cost function is a function D k — > M>o- A valued constraint language is a 
set of cost functions <3? C Func(D, M>o). 

Given a valued constraint language VCSP(<I>) is the problem of taking an instance 
tp, a psm- formula consisting of a sum of m atomic ^-formulas over n free variables x and 
computing the value 

minCost(V') = min /w,(x) , 
xgD" 

where is the function defined by tp. 

We typically use the notation of Kolmogorov and Zivny. An instance is usually denoted 
by the letter /. In this case, we use fj to denote the function specified by the psm-formula 
corresponding to instance /, so the value of the instance is denoted by minCost(I). The 
psm-formula corresponding to I is a sum of atomic formulas (since all of the variables are 
free variables). We refer to each of these atomic formulas as a valued constraint and we 
represent these by the multiset T of all valued constraints in the instance /. For each valued 
constraint t £ T we use ht to denote its arity, ft to denote the function represented by 
the corresponding atomic formula, and at to denote its scope, which is given as a tuple 
(i(t, 1), . . . , i(t, k t )) E {1, . . . , n} kt containing the indices of the variables in the scope. Thus, 

= ^2ft(x i{tA ),...,x i{tM) ). (6) 
teT 

For convenience, we use x[<7t] as an abbreviation for the tuple . . . , x^u^a). In this ab- 

breviated notation, the function defined by instance / may be written //(x) = YLteT ft( x [ a t])- 
Now, let [0, 1]q = [0, 1] n Q. For reasons which will be clear below, it will be useful to 
work with weight functions in Func(Z), [0, 1]q). For such a weight function F, let the cost 
function 1(F) G Func(Z),M>o) be the function defined by 



(*(F))(x) 



-lnF(x) ifF(x)>0 
oo if F(x) = 0. 
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For example, £(EQ) = eq. For a weighted constraint language T C Func(L>, [0, 1]q), let i(F) 
be the valued constraint language defined by 1(F) = {£(F) \ F G F}. 

There is a natural bijection between instances of #CSP(J 7 ) and VCSP(£(J r )), obtained by 
replacing each function Ft in the former by the function ft = (-{Ft) in the latter, keeping the 
scopes unchanged. Note that //(x) = — lni^x), for any assignment x, with the convention 
— InO = oo. 

Definition 17. A valued constraint language is conservative if it contains all arity-1 cost 
functions D — > K>o- 

The mapping F i-> £(F) from Func(D, [0, 1]q) to Func(D, M>o) is not surjective because 
there are real numbers that are not the logarithm of any rational. For the same reason, the 
valued constraint language 1(F) is not conservative (for any weighted constraint language F). 
Finally, note that we have only defined 1(F) for F G Func(D, [0, 1]q). The obvious extension 
to F G Func(D, Q>o) would produce negative-valued cost functions and we wish to avoid this 
since Kolmogorov and Zivny [24] do not allow it. 

Definition 18. A cost function is crisp [13] if /(x) G {0, oo} for all x. 

Definition 19. For any cost function /, let Feas(/) be the relation defined by Feas(/) = {x | 
/(x) < oo}. 

Thus, any cost function / can be associated with its underlying relation. Similarly, we can 
represent any relation by a crisp cost function / for which /(x) = if and only if x is in the 
relation. A crisp constraint language is a set of relations, which we always represent as crisp 
cost functions, not as functions with codomain {0, 1}. For a valued constraint language 3>, 
the crisp constraint language Feas(<3?) is given by Feas(3>) = {Feas(/) [ / G <£}. 

Definition 20. A crisp constraint language is conservative if it includes all arity-1 relations. 

A relational clone is simply a crisp constraint language Feas((<I>}y) for a valued constraint 
language 

Lemma 21. Suppose $ C Func(D, M> ). Then (Feas(<3?))y = Feas((<J>)y). 

Proof. The mapping p: M>o — > {0, oo} defined by p(oo) = oo and p(x) = 0, for all x < oo, is 
a homomorphism of semirings, from (R>o, min, +) to ({0, oo}, min, +). □ 

5 STP/MJN multimorphisms and weak log-supermodularity 

In |24t Corollary 12], Kolmogorov and Zivny give a tractability criterion for conservative 
VCSPs. In particular, they show that the VCSP associated with a conservative valued con- 
straint language $ is tractable iff $ has an STP/MJN multimorphism. 

We define STP/MJN multimorphisms below. In this section, we show (see Theorem 1341 be- 
low) that, if a weighted constraint language T G Func(Z), [0, 1]q) is weakly log-supermodular, 
then the corresponding valued constraint language 1(F) has an STP/MJN multimorphism. 
In Section [6l this will enable us to use such a multimorphism (via the work of Kolmogorov and 
Zivny [24] and Cohen, Cooper and Jeavons [13]) to prove #BIS-easiness and LSM-easiness of 
the weighted counting CSP. 

Our proof of Theorem 1341 relies on work by Kolmogorov and Zivny [23] and Takhanov [29] . 
We start with some general definitions. Most of these are from [23], but some care is required 
since some of the definitions in [23] differ from those in |13j . 



14 



Definition 22. A k-ary operation on D is a function from D to D. An operation on D is 
a fc-ary operation, for some /c. 

We drop the "on D v when the domain D is clear from the context. 

Definition 23. A fc-tuple (pi, . . . , p&) of fc-ary operations pi,...,pk is conservative if, for 
every tuple x = (xi, . . . ,x&) G D fc , the multisets {{xi, . . . ,Xk}} and {{pi(x), . . . , pfc(x)}} are 
equal. 

Note that we have now defined conservative operations and conservative constraint lan- 
guages (weighted, valued and crisp). There are connections between these notions of "con- 
servative" but we do not need these here. 

Definition 24. (pi,...,pfc) is a multimorphism of an arity-r cost function / if, for all 
x 1 , . . . , x fc G D r , we have: 

k k 

i=i i=i 

Definition 25. (pi, . . . , p&) is a multimorphism of a valued constraint language $ if it is a 
multimorphism of every / G 3>. 

These definitions imply the following. 

Observation 26. If (pi, . . . , p^) is conservative, then it is a multimorphism of every unary 
cost function /. 

Definition 27. Suppose M C D 2 . A pair (n, U) of binary operations is a symmetric tourna- 
ment pair (STP) on M if it is conservative and both operations are commutative on M. We 
say that it is an STP if it is an STP on D 2 . 

Definition 28. Suppose M C D 2 . A triple (Mj 1, Mj2,Mn3) of ternary operations is an MJN on 
M if it is conservative and, for all triples (a, b, c) G D 3 with {{a, b, c}} = {{x, x, y}} where x 
and y are distinct and (x, y) G M, we have Mj l(a, b, c) = Mj2(a, b,c) = x and Mn3(a, 6, c) = y. 

Definition 29. An STP/MJN multimorphism of a valued constraint language $ consists 
of a pair of operations (n,U) and a triple of operations (Mj 1, Mj2, Mn3), both of which are 
multimorphisms of for which, for some symmetric subset M of D 2 , (n, U) is an STP on M 
and (Mjl,Mj2,Mn3) is an MJN on {(a, b) G D 2 \M \a^b}. 

Definition 30. $ C Func(D, ]R>o) is weakly submodular if, for all binary functions / G {&)v 
and elements a,b G D, 

f(a,a) + f(b,b)<f(a,b) + f(b,a) or f(a, a) = f(b, b) = oo. (7) 

Note that the definition of weak submodularity for cost functions is a restatement of 
Kolmogorov and Zivny's "Assumption 3". It is not trivial that weak log-supermodularity 
for T is related to weak submodularity for i(J-). Expressibility for VCSP is different from 
expressibility for #CSP and, specifically, we cannot expect (l{F))v = ^((^ r )#) to hold in 
general. However, the following is suitable for our purposes. 
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Lemma 31. Suppose J- C Func(-D, [0, 1]q) and let = £(J~). If T is weakly log-supermodular 
then $ is weakly submodular. 

Proof. We prove the contrapositive. Suppose / G (3>)y is a binary function that witnesses 
the fact that <3? is not weakly submodular according to Definition 1301 specifically, 

f(a,a) + f{b,b) > f(a,b) + f(b,a) and min{/(a, a), f(b, b)} < oo. 

Since / G (3>)y, we may express / in the form 



/(x) =min5f(x,y) = min V ^(x, y), 
y y r— ' 

where the ^ G $ are atomic. For fc G N, define 

m 

FM(x)^n^( x .y)> 

y i=l 

where each Gj is such that gi = £(G{). Note that G (J 7 )#, and 

F (*)( x )iA 

-> max JjGj(x,y), as A; ^ oo. 



y i=i 



Now 



max Gj (x, y ) = max exp ( - ^ ^ (x, y ) J 

i=l ' ^ i=l ' 

= exp ^ - min ^ gi (x, y)^j 



i=l 

exp(-/(x)) 



and 



exp(-/(a, a)) exp(-/(6, 6)) < exp (-/(a, 6)) exp(-/(6, a)). 

Thus F(a,a)F(b,b) < F(a,b)F(b,a) where F = for some sufficiently large k. Also, 
min{/(a, a), 6)} < oo implies that max{F(a, a), F(6, b)} > 0. These properties of F 
imply that T is not weakly log-supermodular, according to Definition SI □ 

Let r be a crisp constraint language. A majority polymorphism of T is a ternary oper- 
ation p such that p(a, a, 6) = /)(a, 6, a) = /)(&, a, a) = a for all a,b £ D and for all arity-/c 
relations i? G T we have 

x,y,z G i? ==> (p(x 1 ,y 1 ,z 1 ),...,p(x k ,y k ,z k )) eR. 

Let N(a,b,c,d) be the relation {(o, c), (6, c), (a, d)}. The existence of such a relation in 
(r)y indicates that V is not "strongly balanced" in the terminology of [20J. Note that, on the 
Boolean domain, N(0, 1,0, 1) is the "NAND" relation. 

Theorem 32. (Takhanov) Let T be a conservative relational clone with domain D. At least 
one of the following holds. 
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• There are distinct a,b £ D such that N(a, b, a, b) € T. 

• There are distinct a,b £ D such that {(a, a, a), (a, b, b), (b, a, b), (6, b, a)} € T. 

• For some k > 1, there are ao, • • • , a 2 k, bo, ... , b 2 k £ -D suc/i i/mf, for each < i < 2k, 
chi 7^ h and, for each < i < 2k — 1, 

JV(ai,&i,Oi + i,&i + i) € T and N(a 2k ,b 2k ,ao,b ) € T. 



• r /ias a majority polymorphism. 

Proof This formulation is essentially [29, Theorem 9.1] except for the last bullet point. As 
stated in the proof of that theorem, the first two conditions both fail if and only if the 
"necessary local conditions" |29l Definition 3.5] hold. Unfortunately for us, Takhanov uses 
the term "functional clone" differently to how we use it, so the reader will need to take 
this into account to understand the local conditions. However, we do not need the detail, 
here. Takhanov's proof of the NP-hard case of his Theorem 3.7 (at the end of his Section 4) 
shows the following: Given the necessary local conditions, the third condition fails only if a 
certain graph Tp is bipartite. If Tp is bipartite then [29, Theorem 5.5] establishes a majority 
polymorphism. □ 

Lemma 33. If $ C Func(D, M>o) is conservative and weakly submodular, T = (Feas(<I>)}y 
has a majority polymorphism. 

Proof. Since is conservative (Definition \T7\i , so is T (Definition [20]) . We will now show that 
the first three bullets of Theorem [32] contradict the premise of the lemma, so the fourth must 
hold. 

The first bullet-point is easily ruled out. Suppose the given relation is in F By Lemma [2TT 
there is a binary function / € {<&)v such that Feas(/) = N(a,b,a,b). This function has 
fib, b) = oo and f(a, a), f(a, b), f(b, a) < oo, and hence violates ([7]). 

For the second bullet-point, by Lemma |2"T1 we have a 3-place function g G (&)v which 
is finite precisely on {(a, a, a), (a, b, b), (b, a, b), (b, b, a)}. Now let M be a sufficiently large 
constant and let u be the unary function defined by 



u(z) 



M \i z = a, 
if z = b, 
oo otherwise. 



Let 



f[x, y) = min{g(x, y, z) + u(z)}. 



Then f(a,a) = M + g(a, a, a), f(b,b) = M + g(b, b,a), f(a,b) = g(a,b,b) and f(b,a) = 
g(b,a,b). Clearly, / violates ([7]) for sufficiently large M. 

Finally, let us consider the third bullet-point. By Lemma [211 we have binary functions 
go, g\, . . . , g 2 k € ($)y where the underlying relation of Oj is iV(oj, bi, aj + i, for < i < 2k 
and the underlying relation of g 2 k is N(a 2 k, b 2 k, a-o, bo). Define 

f(x,y) = mm{g (x,z ) +u (^o) +9i(zo,Zi) + u x (zi) + 

— v u 2k ~i{z 2k -i) + g 2 k(z 2 k-i,y) I (zo,...,z2k-i) e D 2k }, 
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where Uj(aj) = M, itj(&i) = 0, and Ui(z) = oo if z ^ {aj,6j}. Note that /(ao,ao) > fcM, 
/(&0) ^o) > (& + 1)M and /(ao, 60), /(^o, ao) < + (2/c + l)m, where m is the largest finite 
value taken by any of go, . . . ,g%k- So / violates (J7|) for sufficiently large M. 

So we are left with the remaining possibility that T has a majority polymorphism. □ 

We can now prove the main result of this section. 

Theorem 34. Let J 7 be a weighted constraint language such that Funci(Z), [0, 1]q) C J C 
Func(D, [0, 1]q). If T is weakly log-supermodular theni(J-) has an STP/MJN multimorphism. 

Proof. Let $ = £(J")UFunci(D,R> ). We will show that $ has an STP/MJN multimorphism. 
By Definitions 1291 and 1251 this is also an STP/MJN multimorphism of the subset £(J-). 

By Lemma EH 1(F) is weakly submodular. Now, £(7-") contains £(Funci(Z?, [0, 1]q)). Thus, 
for every unary function u G Funci(Z), R>o) and every e G (0,1), there is a unary function 
u £ G £(J 7 ) such that, for all x G {0, 1}, \u(x) — u £ (x)\ < e. From the definition of valued clones, 
and continuity, we deduce that, for every binary function / G (3>)y and every e > (0, 1), there 
is an f £ G (^(J 7 ))^ such that, for all x,y £ {0, 1}, \f(x, y) — f £ (x, y)\ < e. Since i[T) is weakly 
submodular, we conclude from the definition of weak submodularity (Definition I50"j) that $ is 
weakly submodular. 

In [244 §6.1-6.4], Kolmogorov and Zivny show how to construct an STP/MJN multimor- 
phism of $ under "Assumptions 1-3" . Assumption 1 is that $ is conservative, which is true 
by construction. Assumption 3 is that $ is weakly submodular. This is given as a premise of 
our lemma. Assumption 2 is that T = Feas(<I ) ) has a majority polymorphism, which follows 
from Assumptions 1 and 3 by Lemma [331 (Assumption 2 states that $ has a majority poly- 
morphism. In our terminology, this means that Feas($) has a majority polymorphism.) □ 

6 LSM-easiness and #BIS-easiness 

Our aim is to show that if £(J~) has an STP/MJN multimorphism then J- is LSM-easy. This 
will involve using the arguments of [13] and [23], but we try, as much as possible, to avoid going 
into the details of their proofs. We start by generalising the notion of an STP multimorphism. 

Definition 35. Let / be an arity-fc cost function. A generalised STP multimorphism of f 
is a pair (n, U), defined as follows. For 1 < i < k, rij and Uj are operations on the set 
Di = {a £ D \ 3x : Xi = a and /(x) < 00}, and (rij, Uj) is an STP of {/}. 

The operation n is the binary operation on D± X ■ ■ ■ X defined by applying rii, . . . , 
component-wise. Similarly, U is defined by applying Ui, . . . , U/% component- wise. We require 
that, for all x, y G D k , /(U(x, y)) + /(n(x, y)) < /(x) + /(y). Equivalently, we require 

/(□i(xi,yi),. . . ,U k (x k ,y k )) + f(n 1 (x 1 ,y 1 ), . . . , n k (x k , y k )) < /(x) + /(y). 

Where it is clearer, we use infix notation for operations such as n and U. 

Theorem 36 (Kolmogorov and Zivny). Suppose $0 is a finite, valued constraint language 
which has an STP/MJN multimorphism. Then there is a polynomial-time algorithm that 
takes an instance I of VCSP^o) and returns a generalised STP multimorphism (n, U) of 
fj. The pair (n,U) depends only on the STP/MJN multimorphism of an d on the relation 
Feas(/j) underlying fj. It does not depend in any other way on I. 



18 



Proof. This theorem is closely related to Theorem 11 in [23]. To prove Theorem [36], we just 
take the proof of Kolmogorov and Zivny's Theorem 11 (refer to [24, §7]), but we stop before 
the final step. Specifically, we stop when their Lemma 35 has been proved. At this point, the 
existence of the pair (n,U), as specified above, has been established. 

Note that Kolmogorov and Zivny restrict to rationals, whereas we allow real numbers. 
This is not a problem because their algorithm does not require access to the functions in <l>o 
themselves. Instead, it only requires access to the relations in Feas(<^o) an d to the STP/MJN 
multimorphism that <3?o satisfies. These are both finite amounts of data, which can be hard- 
wired into the algorithm, whose input is just the instance /, which is a symbolic expression. □ 

We will also use the following algorithmic consequence of [24| Theorem 1 1] . We restrict to 
crisp cost functions because this is all that we use and we wish to avoid issues with number 
systems. 

Theorem 37 (Kolmogorov and Zivny). Suppose that <3?o *s a finite, crisp constraint lan- 
guage that has an STP/MJN multimorphism. Then there is a polynomial-time algorithm for 
VCSP($ )- 

For our eventual construction, we would like U) to induce a generalised STP mul- 
timorphism of ft for each individual valued constraint t in the instance. We do not know 
whether this is true of the generalised STP multimorphism provided by Kolmogorov and 
Zivny's algorithm, but something sufficiently close to this is true. 



Definition 38. For an instance /, a valued constraint t and a length-fet vector a, define 



if there exists x with x[<t<] = a and //(x) < oo; 
otherwise, 



and define f[ = ft + Rj , 



Thus, f[ is a "trimmed" version of ft whose domain is precisely the ^-tuples of values 
that can actually arise in feasible solutions to instance /. We will see that if the scope at 
contains variables with indices i(t, 1), . . . , i(t, kt), then 

<n[cr t ], u[<7t]} = ((rii(t,i), • • • , r\tM))> ( u i(*,i)' • • • > u i(tM))) 

is a generalised STP multimorphism of //, even though it might not necessarily be a gener- 
alised STP multimorphism of ft- 

Note that Theorem 1371 has the following consequence. 

Corollary 39. Let <3?o be a finite, valued constraint language that has an STP/MJN multi- 
morphism. There is a polynomial-time algorithm that takes an instance I of VCSP(3>o)> a 
valued constraint t and a length-kt vector a and returns a truth table for f[. 

The truth table produced by the algorithm of Corollary [39] is finite since all valued con- 
straints in $0 are finite. 

Theorem 40 (An extension to Theorem [36]) . Suppose $0 is a finite, valued constraint lan- 
guage which has an STP/MJN multimorphism. Consider the algorithm from Theorem [36] 
which takes an instance I of VCSP(<3?o) (in the form ^) and returns a generalised STP 
multimorphism (n,U) of //. Then, for all t £ T, (p[at],\J[(Jt\) is a generalised STP multi- 
morphism of f[. 
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Proof. Focus on a particular valued constraint t of /. Let k = kt be the arity of ft, and for 
brevity denote l~l [at] and U [at] by n' and U', respectively. Without loss of generality assume 
at = (1,2, ... , k). We wish to show that 

//(a n' b) + //(a U' b) < //(a) + //(b) (8) 

for all a, b € D\ x • • • x D^. If either //(a) = oo or //(b) = oo, then we are done. Otherwise, by 
construction of //, there exist x and y such that a = x.[a t ], b = y[at], and //(x), //(y) < oo. 
Notice that //(a) = /t(a) < oo and //(b) = ft(h) < oo, also by construction of //. Now 
consider the augmented instance Jjv of / with N extra copies of the valued constraint t. We 
have 

/ Jjv (x) = /j(x)+iV/ t (a) 

// N (y) = //(y) + iv/t(b) 

/ /iV (x n y) = //(x n y) + Nf t (a n' b) 
/ % (x U y) = / 7 (x U y) + Nf t (a W b). 

Since Feas(// Ar ) = Feas(/j), from Theorem l36l (n, U) is also a generalised STP multimorphism 
of f lN , i.e., 

/ /jv (x n y) + / /jv (x U y) < / Jjv (x) + / 7jv (y). 
Combining this with Q, we obtain 

/ t (a n' b) + / t (a U' b) + OiN- 1 ) < / t (a) + / t (b) + O^ 1 ). (10) 

Since (|10p remains true as iV — > oo but /t is independent of N, we conclude that 

/ t (a n' b) + / t (a U' b) < / t (a) + / t (b). 

Since an'b and aU'b extend to feasible solutions xl~ly and xUy, it follows that //(an'b) = 
ft(a n' b) and //(a U' b) = / t (a U' b). The required inequality © follows immediately. □ 

To make use of Theorem 1401 we wm use the following definitions. 

Definition 41. Given a finite, valued constraint language 3>o C Func(D, M>o), let $g be the 
set of functions of the form / + R, for / € <l>o n Funcfc(Z), R>o), R G Funcfc(Z), {0, oo}) and 
jb g N. 

Note that <I>q is finite because Funcfc(-D, {0, oo}) is finite for any finite k. 

Definition 42. Two n-variable instances I and /' of VCSP( < I , ) are equivalent if //(x) = //'(x) 
for all x € D n 

Definition 43. |13j A function / : D\ x • • • x Z) r — > M>o is domain-reduced if, for each 
i £ {1, r}, and for each a G -Dj, there is an x G -D n such that Xi = a and /(x) < oo. 

Lemma 44. Suppose $0 ^ s a finite, valued constraint language which has an STP/MJN 
multimorphism. Consider an instance I of VCSP^o)- There is an equivalent instance I' of 
VCSP(<3?q) and a generalised STP multimorphism (n, U) of fji which induces a generalised 
STP multimorphism of ft for each valued constraint t of I'. Both I' and (n, U) are polynomial- 
time computable (given I). Moreover, each operation n, and Uj induces a total order. 
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Proof. We first show how to construct an equivalent instance /' and a generalised STP mul- 
timorphism of /// which induces generalised STP multimorphisms on the valued constraints. 
To obtain I', start from the instance / and use Corollary [39] to replace each valued con- 
straint /i(x[<7t]) with //(x[(7i]). This operation clearly preserves the set of feasible solutions 
and their costs. Then use the algorithm from Theorem 00] to construct the generalised STP 
multimorphism (n,U). 

In the remainder of the proof, we construct a new generalised STP multimorphism by 
modifying (n, U) to ensure that it is composed of total orders, as required. Consider the 
following claim. 

Claim: Suppose that D is a domain. Given a set of functions $ C Func(D, R>o), let 
V be an instance of VCSP(<i>) with variable set {v\, . . . , v n }. Let Di = {a € D | 3x : 
Xi = a and /p(x) < oo}. Suppose that (n,U) is a generalised STP multimorphism of V. 
Then there is a generalised STP multimorphism (n' U') ofV in which each is a total order 
on Di (hence is the reversal of this total order). Furthermore, for any set J = . . . C 
{1, . . . , n} and any domain-reduced function f : x ••• x D^. — > M>o for which (n j, U j) 
is a multimorphism, (rij,Uj} is also a multimorphism of f. The multimorphism (rY, U') is 
polynomial-time computable. 

This claim is proved (but not explicitly stated) in the proof of [13j Theorem 8.2]. The 
basic method is as follows. V is augmented with extra (redundant) valued constraints using 
unary and binary crisp cost functions. The binary crisp cost functions are used to enforce 
consistency so that when rij is modified to form a total order on Di, a compatible modification 
is made to each other F\j. Once n and U are constructed, it is proved by induction that every 
relevant function / has the property specified in the claim. The induction is on the arity of /. 

To prove the corollary, we use the claim with f> = $q, V = V and, for each valued 
constraint t of I', f = f[. □ 

Lemma 45. If T C Func(D, [0, 1]q) and has an STP/MJN multimorphism, then 

#CSP((y) is LSM- easy for every finite Q C T . 

Proof. Let Q be a finite subset of T . To any instance I# of #CSP(<?) there corresponds 
an instance / = of VCSP^o); where <I>o = 1(G): for each weighted constraint t, the 

function Ft is mapped to ft = t{Ft) while the scope at remains unchanged. Using Lemma |4"H 
we may construct an equivalent instance I' of VCSP(<3?q) on the domain D\ x • • • x D n and a 
generalised STP multimorphism (n, U) of that instance, where each rij is a total order, (n, U) 
induces a generalised STP multimorphism of each ft- 

We now construct an instance I" over the Boolean domain that is equivalent to I' and 
hence to /. For each i, 1 < i < n, introduce a set of \Di \ + 1 Boolean variables Vi = {zi <a \ a G 
D^}, where Df = Di U {!-}. Extend the total order on Df by placing ± below all elements 
of Di. Define a nested sequence of subsets of Df by Ui ja = {b € Df \ b < a}. The idea is 
that each domain element a G Di is represented by the truth assignment that assigns 1 to all 
variables in Ui a , and to the others. Consider the constraint asserting that only these \Di\ 
particular assignments to Vi are allowed. This constraint can be represented by the crisp cost 
function / that assigns /(x) = to these assignments and /(x) = oo to all others. Note that 
F(x) = exp(— f(x)) is log-supermodular. 

Note that we can use the same relation for any pair of sets Di and Dj with \Di \ = \DA - 
if Di and Dj have different total orders then the relation is applied to the variables in Df in a 
different order than to the variables in Df . If we add these crisp valued constraints then there 
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is a natural bijection between D\ x • • • x D n and feasible assignments to Boolean variables 
V\ U • • • U V n . The variable where a is the smallest (respectively largest) element of Df 
always takes on the value 1 (respectively 0), and so these variables are redundant. However, 
their introduction simplifies the description of some constructions later in the proof. 

Consider a valued constraint in I 1 of arity k that imposes the function /' € ^q, and, 
without loss of generality, assume that its scope is the first k variables x±, . . . ,x k . Add a 
corresponding valued constraint /" to I" with /": 2 lU " k — > M>o defined as follows, where 
for convenience we are viewing /" as a function on subsets of V\ U ■ • • U V k rather than as a 
function of |Vi| + • • • + \V k \ Boolean variables: 



f"(A) 



f'(ai, . . . ,afc), if A = ?7 1)01 U • • • U U kjak for some (ax,. . .,a k ); 
oo, otherwise. 



We claim /" is submodular, i.e., f"(AnB) + f"(AUB) < f"{A) + f"{B). If either f"(A) = oo 
or f"(B) = oo there is nothing to prove. So A = U\ >ax U • • • U U k;Clk and B = Ui tbl U • • • U U k> b k 
for some (a±, . . . , a k ), . . . , 6fe) £ D\ x • • • x D k . Then 

f"(A DB) + f"(A U B) 

= f"((Ui, ai n u 1M ) u • • • u (u ktah n u k , bk )) 

+ f"{(Ui, ai u u lM ) U • • • U (c/ fc , afc u u Kbk )) 

= f"(U lt a inibl U • • • U £/fc,a fe n fc & fc ) 
,aiUibi 

u • • • u c/ fc 

= f'(ai III h, . . . , a k n& + /'(ai Ui 6i, . . . , a k U k b k ) 
< f'(ai,---,a k )+f'(h,...,b k ) 
= f"(A) + f"(B). 

Now take stock. We have an instance I" of Boolean VCSP, which is equivalent to /' and 
hence to /. It has at most n(|D| + 1) Boolean variables and it has n more valued constraints 
than /. The number of distinct valued constraints in $q is | <&q | < | $q | + | D \ ; note that these 
come from a fixed set of cost functions independent of the instance I and hence of Is, itself. 

Now map the VCSP instance I" back to #CSP to yield an instance VL over the Boolean 
domain in which every valued constraint comes from a certain fixed set of cost functions 
J% C LSM. Specifically, I" = £(I^) and $g = £(J%). Since I" is equivalent to I, there is a 
bijection between the non-zero terms of Z(I#) and that preserves weights, and hence 

Z(I # ) = Z(I£). □ 

Lemma 05] shows that, if £(J~) has an STP/MJN multimorphism, then #CSP(£?) is LSM- 
easy for every finite Q C T . Lemma H51 b elow strengthens the result by showing that #CSP(£/) 
is #BIS-easy. The strengthening applies when the weight functions in J- have arity at most 
two. 

In order to do the strengthening, we need to generalise the notion of a binary submodular 
function to cover binary functions over larger domains. A matrix M has the Monge property 
if, for every pair of rows r and r' and every pair of columns s and s', 

M r hr' sAs' + M r y r ' s y s i < M r/ \ r ' s Vs' + M r y r ' s /\ s ' , 

where the A and V operators are the minimum and maximum operators, respectively. To 
apply this concept here, suppose that / is a function with domain Di x Dj. Given orders 
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on Di and Dj, let Di(£) and Dj(£) denote the £'th element of Di and Dj, respectively. We 
say that a function / with domain Di x Dj is Monge (with respect to the given orders) if 
the \Di\ x \Dj\ matrix Mj is Monge), where, as in Section [31 {Mf)m = F(Di(k), Dj(£)). 
We say that a function F : D{ x Dj — )■ [0, 1]q is log-anti-Monge (with respect to the given 
orders) if the function -£(.F) is Monge (with respect to the same orders). The "anti" in the 
term "log-anti-Monge" comes from the fact that (£(F))(x) = — lni^x), and the minus sign 
reverses the inequality in the definition of Monge. Thus, the relationship between Monge and 
log-anti-Monge is analogous to the relationship between submodular and log-supermodular. 
The Monge property generalises submodularity and the property of being log-anti-Monge 
generalises the notion of log-supermodularity. 

Lemma 46. If J- C Func(D, [0, 1]q) is a weighted constraint language whose weight functions 
have arity at most two and £(F) has an STP/MJN multimorphism, then #CSP(<5) is #BIS- 
easy for every finite Q C T . 

Proof. Let Q be a finite subset of T . We use exactly the same construction as in the previous 
lemma, but go further and show that every weight function F" appearing in instance VI 
is expressible in terms of unary weight functions in £Y{ ,i}> an d the binary weight function 
IMP defined by IMP(0,0) = IMP(0, 1) = IMP(1,1) = 1 and IMP(1,0) = 0. Moreover, 
unary weight functions in lA^ x] (even those taking irrational values) can be approximated 
sufficiently closely by polynomial-sized pps- formulas using IMP Lemma 36]. This will 
complete the proof, since #CSP(IMP) <ap #BIS by [17, Theorem 5]. 

The task then, is to show that every weight function F" in instance I'L is expressible 
in terms of unary weight functions in ^{0,1} an d IMP. We do this by considering, in turn, 
the different types of weight functions arising in I'L The n relations (crisp cost functions) 
that were introduced in I" to impose a total order on the variables in the sets Vi are clearly 
implementable in terms of imp = ^(IMP). 

Every other weight function F" is associated with a cost function /" in I" that is an 
implementation over the Boolean domain of a cost function /' from I'. Since /' G $0, it has 
arity at most 2. Our goal is to show that the function -F'(x) = £ (/'(x)) = exp(— /'(x)) 
is expressible in terms of unary weight functions in Z^{o,i} an d IMP. If /' is unary, this is 
immediate, so suppose /' is binary. 

To fix the notation, suppose that /' is a function /': Di x Dj — > M>o- We can assume 
without loss of generality that Di and Dj are disjoint (otherwise, rename some elements). 
Also, Di and Dj are ordered according to the linear order induced by n. Since (n, U) induces 
a generalised STP multimorphism of /' (see Definition l35j) , the function /' is Monge (with 
respect to this order). We start by considering two special cases. 

• First, suppose that F' is strictly positive. (That is, there is no (x, y) with f'(x,y) = 00 
so the range of /' is contained M>o)- Rudolf and Woeginger [28] have shown that every 
nonnegative Monge matrix is expressible as a positive linear combination of certain 
simple Monge matrices. Translated to our setting by applying £~ 1 , this says that F' 
is expressible as a product of certain simple basis functions, namely: (i) the unary 
functions 

Ba ( X ) = h ifx=a; 

I 1, otherwise, 
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for all a G D{ (with a similar set of unary functions defined over Dj), and (ii) the binary 
functions 

a, if x < a and y > b; 
1, otherwise, 



B a ,b{x,y) 



for all (o, b) G (Di, Dj) (with a similar set of binary function defined by replacing x < a 
and y > b by x > a and y <b). In both cases (i) and (ii) , a is an arbitrary constant in 
the range (0, 1). 

Using the Boolean variables Vi U Vj, the basis function B a {x) may be implemented as 
o(^i,a~); where 

V, if 2 = 0; 
1, if 2 = 1, 



Up(z) 



and a is the element immediately below a in the total order on Di. Also, the basis 
function B a , b (x,y) may be implemented as lMP(z ija ,y)TMP(zj ib - ,y)U a (z i)a ) Ui/ a _i(y), 
where y is a new variable. 

• Second, suppose that the range of F' is {0, 1} (that is, /' is a crisp cost function). As 
in the proof of Lemma [431 let U^ a = {b G Df \ b < a}. Let F* be the set of subsets of 
Vi U Vj defined as follows. 

F* = {Ui :Cl U Ujf, I F'(a, b) = 1} 

F' is log-anti-Monge so 

F'(x A x', y A y')^'^ V x', y V y') > F'(x, y)F'{x', y'). 

Thus, F* is closed under unions and intersections. Hence by [18\ Corollary 18], F* is 
a conjunction of implications and constants, and hence can be implemented using IMP 
and unary weights. 

To finish the proof, we use a construction from Lemma 23]. Assume without loss of 
generality that Di = {1, . . . , \Di\} and Dj = {1, . . . , \Dj\} with the usual total order. Define 
H' : Di x Dj — > [0, 1]q by putting H'(x,y) = F'(x,y) for all x G Di, y G Dj except that, if 
F'(l, 1) = 0, we put H'(l, 1) = 1. Note that H' is log-anti-Monge (by construction, since F' 
is log-anti-Monge) and H'(l, 1) / 0. Define G: D- t x Dj -> Q> by 

G(x,y) = max{#'(x* y*) n x+v ~ x * ~ y * \ x* < x and y* < y} , 

where 

/i = max{m > \H\x*, y*)m x+y ~ x *- y * < H'(x, y) 
for all (x,y), (x*,y*) G A x L>j 
with -ff'(x, y) > and x* < x and y* < y}. 

From the definition of fi, we can see that /x > 0, so G is strictly positive. We will show 
that G is log-anti-Monge. Let x, x' G -Dj and y, y' G -Dj. There exist x* < x and x'* < x' and 
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y* < y and y'* < y' such that 

G(x,y)G(x',y') = H'(x*,y*)H'(x'*,y'*){j, x+ y +x ' +v '- x *- y *- x '*- 



y' 




Now, from the definition of G and [i we can see that H'(x,y) = RH>(x,y)G(x,y), where Rh 1 
is the support of H' (considered as a zero-one valued function). Thus, by the definition of H', 
F'(x,y) = Rpi(x,y)G(x,y). Then, by the above constructions, F' can be implemented in 
terms of IMP and unary functions. (The first of the cases above shows how to implement G, 
which is strictly positive; the second shows how to implement R, whose range is {0, 1}). □ 

To use Lemmas 1451 and l36l we need to perform some scaling. For any fc-ary weight function 
in F € F, let mj? = max{/(x) | x G D k }. Let 



and let A(J 7 ) = {A(F) | F € F}. Note that A(F) always takes values in [0, 1]q and that, 
since F is conservative, Funci(L>, [0, 1]q) C A(J r ). 
We return, once more, to the proof of Theorem [6] 

Theorem 47. Let F be a weakly log-supermodular, conservative weighted constraint language 
taking values in Q>o- 

• For any finite Q C T, there is a finite Q 1 C LSM such that #CSP(£) < A p #CSP(a')- 

• If F consists of functions of arity at most two, then #CSP(G) is #BIS-easy for any 
finite Q C T . 

Proof. By Theorem 1341 l(X(F)) has an STP/MJN multimorphism. The result follows from 



Theorem [6j our classification of the complexity of approximating #CSP(J r ), now follows 
from Theorems [TU1 [TBI and E?! 

7 Algorithmic aspects 

Finally, we consider the algorithmic aspects of the classification of Theorem [6j Intuitively, 
there is an algorithm that determines the complexity of #CSP with constraints from a finite 
language % plus unary weights because weak log-modularity is essentially equivalent to bal- 
ance and weak log-supermodularity is essentially equivalent to the existence of a STP/MJN 
multimorphism. Balance and the existence of STP/MJN multimorphisms depend only on 
certain finite parts of the weighted constraint language so balance is decidable by [9] and the 
existence of STP/MJN multimorphisms can be determined by brute force, or by using more 
sophisticated methods from [23]. 

We need to determine whether the infinite language % U Ud is balanced. Fortunately, it 
suffices to check whether 7i UW^ is balanced, where U' D = Funci(.D, {1, 2}), which is finite. 




F/mp if mp > 1 
F otherwise 



Lemmas [35] and [36] and the fact that #CSP(J") < A p #CSP(A(J 7 ))- 



□ 
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(Note that it is not enough to test whether Ti is balanced; also, there is nothing special about 
1 and 2: any pair of distinct, positive rationals would do. In fact, \W D \ = 2\ D \ and there are 
sets of size \D\ which would suffice, but we do not need this here.) 

Lemma 48. Let Ti be a finite, weighted constraint language taking values in Q>o- The 
following are equivalent: (1) H\MA' D is balanced; (2) every finite subset ofHUUn is balanced; 
and (3) HUIAd is balanced. 

Proof. (2) and (3) are equivalent because any pps-formula contains only a finite number of 
atomic formulas. (2) trivially implies (1), since IA' D is finite. It remains to show that (1) 
implies (2) so, towards this goal, suppose that Ti ^W D is balanced. We must show that every 
finite subset ofHUUo is balanced. Suppose that such a subset contains r functions in Ud\H- 

Let {Fi, . . . , F r } be unary functions such that Fi(d) = (i 6 {1, . . . , r}, d G D) and 
let Q = % U {Fi, . . . , F r }. We may consider the a^d as formal variables and treat a function 
G £ {G)# with free variables x as a function of both x and the a^- We will show that, for any 
function G and any interpretation of the am (i.e., any instantiation of the function symbols 
Fi as concrete functions D — > Q>o), the matrices associated with G have block-rank 1, thus 
establishing that Q is balanced. 

So, consider any G S (G)# with arity n > 2 and choose any k with 1 < k < n. We will 
show that the D k x D n ~ k matrix Mc(x,y) has block-rank 1 for any value of the a^. By 
Lemma [11] part ([2]), it suffices to show that every 2x2 submatrix induced by rows x, x' and 
columns y,y' has block-rank 1. By Lemma [TT1 part ([1]), this happens if, and only if, every 
such submatrix has rank 1 or at least two zero entries, which happens if, and only if, the 
multivariate polynomial 

p = G(x, y / )G(x / , y)G(x', y') [G(x, y)G(x', y') - G(x', y)G(x, y')] 

is zero for all values of x, x', y, y' and for all values of the a^- (Note that, if the submatrix 
defined by a pair of rows and columns does not have block-rank 1 but has exactly one zero, 
then only one of the four possible choices for x, x', y,y' will make p non-zero.) 

We now fix x,x', y,y' and consider p as a function of just the a^. Our goal is to show 
that (for every choice of x, x', y,y'), p is identically 0. 

Consider first the case where every aj jC ; is a power of two. Here, every atomic formula 
Fi(z) defines the same function as some product U\{z) ■ ■ ■ U^{z) of atomic formulas from IA' D 
so G is equivalent to some function in (HlMA' D )#. But Ti UU' D is balanced by assumption, so 
p = whenever every is a power of two. Therefore, p = over a space that is a product 
of infinite sets. It follows from the Schwartz-Zippel lemma or from [T] Theorem 1.2] that the 
only polynomial with this property is the zero polynomial, so p is the zero polynomial and 
H U {F±, . . . , F r } is balanced for any set {Fi, . . . , F r } of unary weights. □ 

Theorem 49. There is an algorithm that, given a finite, weighted constraint language T~L 
taking values in Q>o, correctly makes one of the following deductions, where T = % \JUd'. 

1. #CSP(t?) is in FP for every finite Q C T; 

2. #CSP(<5) is LSM -easy for every finite Q C T and ^JilS-hard for some such Q; 

3. #CSP(<5) is #BlS-easy for every finite Q C T and #BIS- equivalent for some such Q; 
4- #CSP(<5) is ^SAT-easy for every finite Q C T and ^SAT- equivalent for some such Q. 
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If every function in % has arity at most 2, the output is not deduction [£[ 

Proof. We reduce the problem to determining whether T~L U U' D is balanced, whether £{H) 
has an STP/MJN multimorphism and whether H contains only functions of arity at most 2. 
Balance of finite languages is decidable [9j. An STP/MJN multimorphism consists of two 
operations D 2 — > D and three operations D s — > D, which must have certain easily checked 
properties with respect to each of the functions in £(1-1). Thus, we can determine the exis- 
tence of an STP/MJN multimorphism by brute force, checking each possible collection of five 
operations, or by using the methods of Kolmogorov and Zivny |24J. It is clearly decidable 
whether H contains a function of arity greater than 2. 

By Lemma H5J if H UW^ is balanced, then so is any finite Q C % UUd- Therefore, by 
Lemma [TH #CSP(£?) can be solved exactly in FP so we output deduction[TJ From this point, 
we assume that H U U' D is not balanced. 

Since "KVMA' D is not balanced, nor is HUUd (Lemma l4*8|) . Therefore, HUUd is not weakly 
log-modular (Lemma fl"3|) so there is a finite Q C % UUd such that #CSP(<?) is #BIS-hard 
(Theorem [10]) . 

t(A(H UW D )) has an STP-MJN multimorphism if, and only if, £(A(H)) does (Observa- 
tion [26]), and i(A(H)) is a finite language so we can determine whether it has an STP-MJN 
multimorphism by exhaustive search. If £(A(T-L U Ud)) has an STP-MJN multimorphism, 
then, for all finite Q C A(HUUd), #CSP(<5) is LSM-easy (Lemma H5]) . Since any function in 
A(7~L UUd) is a scalar multiple of some function in % UUd, #CSP(<5) is also LSM-easy for all 
finite Q C HUUd- We output deduction [21 unless every function in % has arity at most 2, 
in which case #CSP(<7) is #BIS-easy for all finite Q C H UUd (Lemma 06]) and we output 
deduction 

On the other hand, if l(A(HUU D )) has no STP-MJN multimorphism, then A(% UU D ) 
is not weakly log-supermodular (Theorem [M]) . Because A is just a rescaling, H UUd is also 
not weakly log-supermodular. Therefore, there is a finite Q C % UUd such that #CSP((/) is 
# SAT-equivalent (Theorem [TOl again) . We output deduction [H □ 
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